TELKOMNIKA Telecommunication, Computing, Electronics and Control 

Vol. 18, No. 2, April 2020, pp. 1072~1078 

ISSN: 1693-6930, accredited First Grade by Kemenristekdikti, Decree No: 21/E/KPT/2018 

DOI: 10.12928/TELKOMNIKA.v1812.15857 O 1072 


Physical security with power beacon assisted 
in half-duplex relaying networks 
over Rayleigh fading channel: performance analysis 


Phu Tran Tin', Duy-Hung Ha’, Luu Gia Thien’, Tran Thanh Trang’ 
'Faculty of Electronics Technology, Industrial University of Ho Chi Minh City, Vietnam 
*Wireless Communications Research Group, Faculty of Electrical and Electronics Engineering, 
Ton Duc Thang University, Vietnam 
3Posts and Telecommunications Institute of Technology Ho Chi Minh City Campus, Vietnam 
National Key Laboratory of Digital Control and System Engineering, Vietnam 


Article Info ABSTRACT 

Article history: In this research, we proposed and investigated physical security with power 
beacon assisted in half-duplex relaying networks over a Rayleigh fading 
Received Jul 24, 2019 channel. In this model, the source (S) node communicates with 
Revised Nov 19, 2019 the destination (D) node via the helping of the intermediate relay (R) node. 
Accepted Dec 4, 2019 The D and R nodes harvest energy from the power beacon (PB) node 
tS in the presence of a passive eavesdropper (E) node. Then we derived 
Keywords: the integral form of the system outage probability (OP) and closed form 
of the intercept probability (IP). The correctness of the analytical of the OP 
Energy harvesting (EH) and IP is verified by the Monte Carlo simulation. The influence of the main 
Half-duplex (HD) system parameters on the OP and IP also is investigated. The research results 

Intercept probability (SP) indicated that the analytical results are the same as the simulation ones. 


Monte carlo simulation 
Relaying network This is an open access article under the CC BY-SA license. 


Corresponding Author: 


Duy-Hung Ha, 

Wireless Communications Research Group 

Faculty of Electrical and Electronics Engineering 

Ton Duc Thang University, Ho Chi Minh City, Vietnam. 
Email: haduyhung@tdtu.edu.vn 


1. INTRODUCTION 

Radiofrequency (RF) energy harvesting (EH) in the wireless communication network is become a 
novel solution for the communication network with battery-limited devices and has attracted massive 
attention in research and industrial directions. The communication network with the battery-limited devices 
or devices and wireless sensors, which are working in the dangerous conditions is the main inside human 
bodies object of RF EH wireless communication network application. This solution can be considered as the 
main one because of carrying both energy and information of the EF, to help the battery-limited devices to 
harvest energy for information transmission to the destination. This technic is called simultaneous wireless 
information and power transfer (SWIPT) [1-11]. Nowadays, the physical layer security (PLS) in EH 
communication cooperative relaying network is popularly studied with considerable interest. The first 
concept of PLS was proposed by authors in [12, 13]. In this paper, the author proposed the secret 
communication between the source and destination nodes with the presence of the eavesdropper channel. 
Furthermore, the cdperative jammer is used for secure the cooperative relaying network by degrading 
the eavesdropper’s channel is proposed and studied in [14, 15]. In this cooperative network, the jammer not 
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only is used to degrade the eavesdropper’s channel but also is helpful for increasing the EH process 
of the energy receiver in the cooperative relaying network. Moreover, a harvest-and-jam (HJ) protocol with 
multi-relay and multi-node in the cooperative relaying network was proposed in [16] to improve the secrecy 
rate of the energy harvesting and information transmission. Also, different secure relay beam-forming 
algorithms for SWIPT were discussed in [17-20]. 

In this research, we proposed and investigated physical security with power beacon assisted in 
half-duplex relaying networks over a Rayleigh fading channel. In this model, the source (S) node 
communicates with the destination (D) node via the helping of the intermediate relay (R) node. The D and R 
nodes harvest energy from the power beacon (PB) node in the presence of a passive Eavesdropper (E) node. 
Then we derived the integral form of the system outage probability (OP) and closed form of the intercept 
probability (IP). The correctness of the analytical of the OP and IP is verified by the Monte Carlo simulation. 
The influence of the main system parameters on the OP and IP also is investigated. The research results 
indicated that the analytical results are the same as the simulation ones. 


2. SYSTEM MODEL 

In Figure 1, the source (S) node communicates with the destination (D) node via the helping of 
the intermediate relay (R) node. The D and R nodes harvest energy from the power beacon (PB) node in 
the presence of a passive Eavesdropper (E) node. Figure 2 draws the energy harvesting (EH) and information 
processing (IT) of the model system. In this protocol, the transmission is divided into blocks of length T, 
which consists of three-time slots. In the first time slot aT (a is the time switching factor, 0<a<1), the S 
and R harvest energy from the PB node. In the remaining intervals time (1-a)T/2, the source S transfers 
the information to R, and R transfers information to D node [21-25]. 
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Figure 1. System model Figure 2. Time switching protocol 


2.1. Energy harvesting phase 
In the first phase, the power beacon will supply the energy for both S and R nodes. 
Hence, the harvested energy at the source and relay can be given as, respectively 


E; = nPpat |hgs|* (1) 
ER = nPgaT |her|* (2) 


where 0 < 7 < 1 is energy conversion efficiency, Pg is the average transmitted power at the power beacon, 
and hgs, hsr are the channel gain of B-S link, B-R link, respectively. The average transmitted power at 
the source and relay nodes can be obtained from (1) and (2), respectively 





= Es _ nPgat|hgs|* _ 2 
S= a=ar/2 G-r KPs|hss| (3) 
Pr = KPglhgrl? (4) 
2na 


where = —. 
1-a 


2.2. Information transmission phase 
In the second phase, the received signal at the relay can be rewritten as 


Vr = hspXs + N, (5) 
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where hsr is the channel gain of S-R link, Xs is the transmitted signal from source and n, is additive white 
Gaussian noise (AWGN) with variance No. In the third phase, the received signal at the destination can be 
given by 

Ya = NppXy + Ng (6) 


where hro is the channel gain of R-D link, x; is the transmitted signal from relay and ng is (AWGN) with 
variance No. Here, we consider amplify and forward (AF) mode at R. Hence, the amplifying factor can be 
given as 


_ Xr PR 
Pim ag 2 (7) 
Yr Pslhsri“+No 


substituting (7) into (6), we have 


Ya = hrpXYyr = MrpxlhspxXs + nr] + Ng = hsrhrDXXs + hrpX%nr + Na (8) 
signal noise 


3. SYSTEM PERFORMANCE ANALYSIS 
From (8), the end to end signal to noise ratio (SNR) of S-R-D link can be calculated as (9). 


y _ Ef{|signal|*} _ |hsr| lhrpl?x?Ps _ lhsrl\lhrpl Ps (9) 
PRD E{|noise|?} lhRDl? X? No +No Ine I7Not+ 2 


After doing some algebra and using the fact that No<<Pr, the (9) can be rewritten as (10). 


= lnsri*lARDI*PsPR 10 
VRD = Sec Pee pone (10) 
Prilhrpl“NotNoPslhsrl 


Substituting (3) and (4) into (10), the SNR can be reformulated as (11). 


_ lAsrl*larpI*«¥ lhpsl*|hpri* _ KYXY 


VSRD = 5-2 (11) 


lhprl*larp|7+lABs|7lAsri7 X+Y 


Where =2, x= lhgrl llron Y = lhgsl |hsr|? . The received signal at the eavesdropper can be 
(0) 


given by 
Ve = lpgeXr + Ng (12) 


where hre is the channel gain of R-E link and ng is AWGN with variance No. The SNR at the eavesdropper 
can be expressed as 


lhpge|*P 
Yp = a (13) 


substituting (4) into (13), we have: 
Ye = KY heel lapel? = KPZ (14) 
Where Z = Kee Feel" 


Lemmal. Please note that all channel are the Rayleigh fading channels, so the probability density function 
(PDF) of |h,;|7 can be given by: 


fing 2) = aye“ (15) 


where i E€ (SR, RD, BS, BR, RE). Moreover, the cumulative distribution function (CDF) of |h;|? also can be 
obtained by 


Fpp) = 1- e74 (16) 
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where/,is the mean value of the exponential random variable |/,|7. 
Lemma 2. The CDF of Xand Y can be computed as respectively 


F(a) = fo Flagyl? (5 lihar? = x) Fjnpg2 dx (17) 


lhBri? 
Fy(b) = fo” Fingal? (a5 Isl? = y) fings EY (18) 
Y o ` lAsrl” \jagsl2 th BPS lABs| 
utilizing the result in [26], the CDF of X and Y can be shown as the below equation, respectively 
Fy(a) = |= 2./AgpAgraKk, (2 årnågra) (19) 
Fy(b) = 1 — 24 AspAgsbky (2y Àsråssb) (20) 


where K„(e°)is the modified Bessel function of the second kind and v" order. From (21) and (22), the PDF 
of X and Y can be calculated as, respectively after applying the formula 


OKn(Z) n 
wen) = Ky (Z) — = Kn (2) 


fx(a) = ee = 2J AppApaho 20 Apnea) (21) 
f(b) = XO = 2,/AspdpsKo(2/AseAnsb) (22) 


3.1. Outage probability (OP) 
The OP of system can be given as 


OP = Pr(Ysrp < Yen) (23) 


where y;,, = 27" — 1 is the threshold of system and R is source rate. Substituting (11) into (23), we have 


OP = Pr (ZŒ < yu) = Pr(XIKYY — yul < Yn) 


Pr (x < -2 ) y > vith 





= KPYY—Yth? ’ KY 
1,Y < 
” T kY 
= fr fy(y)dy + fres F (Gry IY =y) x fr ay (24) 


Applying (19-22), (24) can be rewritten as (25). 





YthY 
ÀsrÀBsÀRDÀBR | x 
OP = i 2./AspAgsKo(2/AsrdAgsy )dy — 4 Sra (25) 
* Ko (2/Asedasy) Ks (2 [Arran [22 ]) ay 





By changing the variable for the first term of (25): t = 2,/AspAgsy. 





Asp ApsArv Apr | rey | x 


KPY-Yth 








1 o0 o0 
OP = ———/f tx K(t)dt — 4 fru, (26) 
VAsRABs `O 3 Ky 
j Ko (24 Asrågsy )Kı (2 ÀRDÅBR | J) dy 
\ th 
Apply (6.561, 16) of the table of integral, (26) can be reformulated as 
thy 
i ÀsrÀBsÀRrDÀBR Fs —| x 
OP = ——.- 4 fy, (27) 





Visries* 
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3.2. Intercept probability (IP) 
The IP can be defined by 


IP = Pr(Vg = Yin) (28) 


substituting (14) into (28) and then applying formulas in (19) or (20), finally we have: 
ApRA ApRA 
IP = Pr(viPZ > ya) = 1 — Pr (Z < te) = 2 [Bestar x x, (2 fPomiaer) (29) 


4. NUMERICAL RESULTS AND DISCUSSION 

The influence of n on the OP and IP is plotted in Figure 3 and Figure 4 for various values of y. Here 
we set a = 0.5, R = 0.5 bps/Hz and y = 1, 5, 10 dB. The first observation one can see from these results is 
that the analytical match very well with the simulation results. Figure 5 and Figure 6 illustrated the optimal 
switching time system OP and IP versus y for n = 0.8, R=0.5 bps/Hz, and a = 0.25, 0.5, 0.85 respectively. 
It should be pointed out that the simulation and analytical results are the same. 


Intercept probability versus 7) with a=0.5 and R=0.5bps/Hz 
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Intercept probability versus Y with 7=0.8,R 
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Figure 6. OP versus y 


5. CONCLUSION 

In this research, we proposed and investigated physical security with power beacon assisted in half- 
duplex relaying networks over a Rayleigh fading channel. In this model, the source (S) node communicates 
with the destination (D) node via the helping of the intermediate relay (R) node. The D and R nodes harvest 
energy from the power beacon (PB) node in the presence of a passive Eavesdropper (E) node. Then we 
derived the integral form of the system outage probability (OP) and closed form of the intercept probability 
(IP). The correctness of the analytical of the OP and IP is verified by the Monte Carlo simulation. 
The influence of the main system parameters on the OP and IP also is investigated. The research results 
indicated that the analytical results are the same as the simulation ones. 
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